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Abstract 

We have performed the first comprehensive global QCD analysis includ- 
ing the CCFR and NuTeV di-muon data; this data provides strong con- 
straints on the strange quark PDF. 

Introduction 

The recent measurements of both neutrino and anti-neutrino production of 
charm (as di-muon final states) by the CCFR and NuTeV collaborations pro- 
vide important new information on the strange quark distribution, s(x), of the 
nucleon.[l,2] We report here the first comprehensive global QCD analysis that 
includes the CCFR and NuTeV di-muon data. 

In previous global analyses, the predominant information on s came from 
the difference of (large) inclusive cross sections for neutral and charged current 
DIS; hence, the comparably small s and s distributions extracted had large 
uncertainties. Lacking better information, these studies often assumed the dis- 
tribution was of the form s{x) = s{x) ^ k{u -\- d)/2 with k, ^ 0.5. The recent 

*The work presented here is performed in collaboration with S. Kretzer, P. Nadolsky, 
W.K. Tung J.F. Owens, S. Kuhlmann, J. Pumplin, D. Stump, M.H. Reno, J. Huston, & 
H. L. Lai. 

t To appear in the proceedings of 11th International Workshop on Deep Inelastic Scatter- 
ing (DIS 2003), St. Petersburg, Russia, 23-27 Apr 2003. 
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Figure 1: a) The basic process v^N ji ji^X . b) The strange quark PDF 
vs. X for the 4 PDF sets discussed in the text. 

high-statistics measurements of cr^^^ and cr^^'^ by the CCFR and NuTeV 
experiments ahow us to separately determine s{x) and s{x) with unprecedented 
accuracy. Neutrino induced di-muon production, v/vN ^^^^X, proceeds 
primarily through the subprocess s — > c or s — > c, and hence provides informa- 
tion on s and s directly. (Cf., Fig la.) 

We present a global analysis including this new di-muon data, corrected 
for experimental cuts and efficiencies using information provided by the exper- 
imental group. The new results are rich in physical content, in part, due to 
the interplay of these high statistics measurements and the strong constraints 
of the PQCD framework. 

Global Analysis 

The inclusion of the CCFR-NuTeV neutrino and anti-neutrino di-muon produc- 
tion data give direct handles on s and s in the important x range ^ [0.01, 0.2]. 
CCFR has recorded 5030 ly and 1060 D di-muon events, while NuTeV has 
recorded 5012 ly and 1458 i> di-muon events. Additionally, NuTeV had a sign- 
selected beam to separate the f and D events. The global QCD analysis contains 
the full data set from the CTEQ6 analysis, in addition to the di-muon data; 
as such, this can be considered as an extension of the on-going CTEQ series of 
global analysis. [3] 

Incorporating this into the global QCD analysis is a non-trivial task. The 
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experimental measurement is the cross section for producing two final state 
muons (dcr^i^T ) with experimental cuts, whereas the theoretical quantities that 
are most directly related to the parton distribution analysis are the underlying 
"charm quark production cross sections" {dacn±) for the process Vn/vjiN — > 
fjr^cX. We can relate these quantities via the equation: 



Here, T and O denote the fragmentation and decay kinematic variables of the 
charmed quark and charmed hadron, respectively; Dc{T) is the fragmentation 
function, and Ac{Q) is the decay distribution function. 

The gap between these two cross sections are usually bridged by Monte 
Carlo programs which incorporate experimental cuts and efficiencies as well as 
fragmentation models. In our analysis, we rely on a Pythia program provided 
by the CCFR-NuTcV collaboration.^ This Monte Carlo calculation in done in 
the spirit and the framework of leading-order (LO) QCD. Accordingly, the the- 
oretical formulas used are also in LO. Since charm production has such a small 
cross section in the experimental kinematic range covered, this approximation 
is perfectly adequate (as we will verify later) for a first study of the di-muon 
data within the global QCD analysis framework. Needless to say, all fully in- 
clusive (large) cross sections used in this study are treated in NLO QCD, as in 
all modern global analyses. 

Results 

We parameterize the strange quark in the form: 



In the present analysis, we take s{x) = s{x); in later studies we shall consider 
a non-symmetric strange distribution. [4-7] 

We will consider four separate sets of PDF's: the CTEQ6M set (which is 
not fit to the di-muon data), a Constrained fit where only the normalization 

^We thank Tim Bolton and Max Goncharov, in particular, for providing this program, as 
well as assistance to use it. Both are vital for carrying out this project. 




(1) 



s{x, Qo) = s(x, Qo) = ao x"' (1 - a;)"=e"^^ (1 + e^^x) 



(2) 
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Figure 2: Plot of {Data — Theory) /Theory for the NuTeV v and D data vs. 
data point number. Note the larger error bars on the D data due to the lower 
statistics. 



{ag} is allowed to vary; a Mixed fit where both {oq, are allowed to vary; a 
Free fit where all {ai} are allowed to vary; 

The CTEQ6M PDF's yield a respectable of 2173 for 1991 data points 
(CTEQ6M was not fit with the di-muons), and the Constrained, Mixed, and 
Free fits yield 2144, 2142, and 2133, respectively. While it is reassuring to see 
the decrease as we increase the number of free parameters, it will require 
additional study to determine whether the data can actually constrain all these 
parameters. 

It is also interesting to note that as we move from the CTEQ6M fit to the 
fits including the di-muons, the x? /DOF for all the experiments (excluding the 
di-muons) are generally unchanged to within a percent. From this observation, 
it is evident that the effect of adding in the di-muon data serves to adjust the 
strange PDF, but has virtually no effect on the other data sets. Or, equivalently, 
this demonstrates that the di-muon data plays a dominant role in determining 
the strange distribution. 

The observation that it is primarily the di-muon data which influence the 
strange distribution prompts us to fit the neutrino and anti-neutrino data sep- 
arately to investigate the extent to which this data can determine the s{x) 
and s{x) distributions independently.^ At leading-order, the neutrino induced 
process vs — > /^^c is tied to the strange distribution, while the anti- neutrino 
induced process vs — > /i^c is tied to the anti-strange distribution. Of course, 

■^A complete analysis which fits the neutrino and anti-neutrino data simultaneously and 
allows six 7^ s{x) is in progress. [6, 7] 
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these relations will be complicated at NLO; however it is useful to see if, and 
how, these separate data sets pull the fit. 

In Fig. lb we display the obtained strange distribution for the combined v 
and V fit, as well as the individual fits. We observe that the v data yields a 
smaller PDF, while the v yields a slightly larger fit than the combined set. Note 
that this effect is consistent with the results of the NuTcV fit by Goncharov,[l] 
which obtained k = 0.35 for the v data and h = 0.41 for the v data. 

This result suggests that the di-muon data is capable of providing informa- 
tion about the a and s distributions separately. However, to properly extract 
this information requires a simultaneous fit to the v and v while imposing the 
strangeness sum rule on the PDFs: / dx s[x) — s{x) = 0. 

Conclusion 

We present the first global analysis which includes the CCFR and NuTeV 
neutrino and anti-neutrino di-muon production data. We find several classes 
of solutions in the strangeness sector that are consistent with all relevant world 
data used in these global analyses. 

While the CTEQ6M PDF provided a good description of the data, we 
obtained an improved fit if we free the parameters of the strange quark PDF to 
allow them to conform to the di-muon data. A more comprehensive study of 
this data is in progress, [6, 7] and future work includes relaxing the constraint of 
s = s, and including the NLO modeling of the charm production cross section. 

We thank members of CCFR and NuTeV collaboration for discussions about 
their experimental results. In particular, we thank T. Bolton for assistance with 
the di-muon Fortran programs. F.O. acknowledges the hospitality of MSU and 
BNL where a portion of this work was performed. This research was supported 
by the U.S. Department of Energy (Contract No. DE-FG03-95ER40908), and 
by the Lightner-Sams Foundation. 
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